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Abstract. We generalize work of Oh & Zumbrun and Serre on spectral stability of spatially periodic traveUng waves 
of systems of viscous conservation laws from the one-dimensional to the multi-dimensional setting. Specifically, we 
extend to multi-dimensions the connection observed by Serre between the linearized dispersion relation near zero fre- 
quency of the linearized equations about the wave and the homogenized system obtained by slow modulation (WKB) 
approximation. This may be regarded as partial justification of the WKB expansion; an immediate consequence is 
that hyperbolicity of the multi-dimensional homogenized system is a necessary condition for stability of the wave. As 
pointed out by Oh & Zumbrun in one dimension, description of the low-frequency dispersion relation is also a first 
step in the determination of time-asymptotic behavior. 



1 Introduction 

Nonclassical viscous conservation laws arising in multiphase fluid and solid mechanics exhibit a 
rich variety of traveling wave phenomena, including homoclinic (pulse-type) and periodic solutions 
along with the standard heteroclinic (shock, or front-type) solutions. Here, we investigate stability 
of periodic traveling waves: specifically, the spectrum of the linearized operator about the wave. 
Our main result generalizes the works |OZ.lj . |Se.lj about stability of periodic traveling waves of 
systems of viscous conservation laws from the one-dimensional to the multi-dimensional setting. 

Consider a system of conservation laws 

ut + Y.f{n),^=Y.{B^'{u)u.,).^, (1.1) 

j 3, k 

u G Z^(open) G M", f G M", B^'' G M"''", eR'^,d> 2, and a periodic traveling wave solution 

u = u{x ■ V — st), (1-2) 
of period X, satisfying the traveling wave ordinary differential equation 

iY^,.,,.kB^\u)uy = i^i^jfiu))' - su', (1.3) 



with initial conditions 
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n(0) = u{X) =: uo- 



Integrating we reduce to a first-order profile equation 

J2 i^jykB^\u)u' = i^jfi^) -su-q (1.4) 

encoding the conservative structure of the equations, where g is a constant of motion. 

The one dimensional study was carried out by Oh & Zumbrun |OZ.lj in the "quasi-hamiltonian" 
case that the traveling- wave equation possesses an integral of motion, and in the general case by 
Serre |Se.lj . An important contribution of Serre was to point out a larger connection between the 
linearized dispersion relation (the function \{^) relating spectra to wave number of the linearized 
operator about the wave) near zero and the homogenized system obtained by slow modulation 
approximation, from which the various stability results of |UZ.lj . |Se.lj may then be deduced. 
The purpose of this paper is to extend to multiple dimensions this important observation of Serre, 
relating the linearized dispersion relation near zero to a multi-dimensional version of the homog- 
enized system developed in |Se.lj . As an immediate corollary, similarly as in |OZ.lj . |Se.lj in the 
one-dimensional case, this yields as a necessary condition for multi-dimensional stability the hy- 
perbolicity of the multi-dimensional homogenized system. In case of stability (so far not found), 
this relation is also the first step in the derivation of asymptotic behavior, as in |()Z.2j in the 
one-dimensional case; this we defer to a future investigation. 

We here make only generic assumptions like those in |Se.lj , ensuring that the set of periodic traveling 
waves is a manifold of maximal dimension subject to the conservative properties of the equations 
(encoded in form H1.4|) '). Given (a, s, u^q) £U X M X S'^-^ X M", ((131) admits a unique local solution 
u{y;a, s^v^q) such that «(0; a, s, i/, g) = a. Denote X the period, u; := 1/X the frequency and M 
and the averages over the period: 

M := ^ r u{y)dy, P := ^ (/^ (n) -Y,B^\u)LOUkdyu)dy 
Jo Jo f^^^ 

when u is a periodic solution of (^3]). Since these quantities are translation invariant, we consider 
the set P of periodic functions u that are solutions of ()1.4() for some triple (s,i^, g), and construct 
the quotient set V := P/TZ under the relation 

(n TZ v) {3h G M; v = u{- - h)). 

We thus have class functions: 

X = X{u), Lo = n{u), s = S{u), V = N{u), q = Q{u), M = M{u), = F^{u), 

where u is the equivalence class of translates of different periodic functions. Note that tt is a 
nonconstant periodic solution. Without loss of generality, assume S{u) = and N{u) = ei, so that 
(|1.4I) takes the form 

B'\u)u' = f{u)-q 
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for q = Q{u). Letting X = X{u) and a = u(0) = uq, the map (y, a, s, u, q) i— s- u{y; a, s, u,q) — a is 
smooth and weh-defined in a neighborhood of {X;a,0,ei,q), and it vanishes at this special point. 
Here and elsewhere, ej denotes the jth standard Euclidean basis element. We assume: 

(HO) P, B^^ G C2. 

(HI) Re aiY^jk ^j^kB^'") > > 0. 

(H2) The map H : RxU xRx S"^-i x M" ^ taking {X; a, s, u, q) I— > u{X; a, s, i/, g) — a is a 
submersion at point (X; a, 0, ei, g). 

As a consequence of (HO), (H2), there is a smooth n+d dimensional manifold V of periodic solutions 
ii in the vicinity of where d is the spatial dimension. On this set, one may obtain, rescaling by 
(x, t) (ex, et) and carrying out a formal WKB expansion as e — > a closed system oi n + d 
averaged, or homogenized, equations 



5tM(u) + ^5,.^,(F^(n)) = 0, 



dt{QN{u)) + \/^{QS{u)) = ^ (1.5) 

in the (n + d)-dimensional unknown u, expected to correspond to large time-space behavior. For 
details, see Section 0] The problem of stability of u may heuristically be expected to be related to 
the linearized equations of H1.5|) about the constant solution u{x,t) = u^, ~ u, provided that 
the WKB expansion is justifiable by stability considerations. This leads to the homogeneous degree 
n + d linearized dispersion relation 

Z^(C,A):=det(A ^ \ u) + Y.'^i du ^^-^^ 

On the other hand, one may also pursue the direct course of linearizing PDE Hl.lf) about the 
stationary solution u and studying the spectrum of the associated linearized operator L. Taking 
the Fourier transform in constant directions Xj, j ^ 1, and following the general construction of 
[U] . |OZ.lj . we obtain an Evans function D{^,X), ^ G M*^, A G C, of which the zero set ('^, A(^)) 
determines the linearized dispersion relation for (jl.lj) . with A(^) running over the spectrum of L 
as ^ runs over W^. For details, see Section [21 In particular, the low-frequency expansion of A(^) 
near (^, A) = (0, 0) may be expected to determine long-time asymptotic behavior, provided that 
spectrum away from A = has strictly negative real part, and this in turn may be expected to derive 
from the lowest order terms of the Taylor expansion of D. A tedious, but fairly straightforward 
calculation following |OZ.lj . |Se.lj shows that 

Z)(e,A) = Z\i(C,A) + 0(|e,Ar2), (1.7) 

where Ai is a homogeneous degree n + 1 polynomial expressed as the determinant of a rather 
complicated 2n x 2n matrix in (^, A): in particular, not in the simple form det AA^o + ^^j^j of 
a first-order hyperbolic (n + 1) x (n + 1) dispersion relation, or an obvious tensor product thereof. 
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Our main result is the following theorem relating these two expansions, generalizing the result of 
|Se.lj in the one-dimensional case. Define 

Z\(e,A):=Ai-^i(e,A), (1.8) 

where A is defined as in (|1.6() . 

Theorem 1. Under assumptions (H0)-(H2), Ai = FqA, i.e., 

Z)(C,A) = roZ\(C,A) + 0(|e,Ar+2) (1.9) 
Iq 7^ constant, for A| sufficiently small. 

That is, up to an additional factor of A"'"^, the dispersion relation H1.6|) for the averaged system 
H1.5|l indeed describes the low-frequency limit of the exact linearized dispersion relation 

^(e,A) = 0. 

The discrepancy A'^"^ is an interesting and at first glance puzzling new phenomenon in the multi- 
dimensional case. However, it is easily explained by a closer look at the formal approximation 
procedure described in Section |^ For, in the derivation of (|1.5j) . it was assumed that QN represent 
the gradient V x<P of a certain phase function (j){x,t). In one dimension, this is no restriction, since 
we may always take (l){x,t) := Lo{z)dz. However, in multidimensions, it imposes the additional 
constraint 

curl (nN) = 0, (1.10) 
which properly should be adjoined to the averaged system. 
Taking the curl of the second equation of H1.5() . we obtain the simple equation 

dt curl (QN) = 0, 

revealing at once that constraint H1.1U() is compatible with the time-evolution of the system, and that 
the unconstrained system possesses {d — 1) spurious zero characteristices A(^) = 0, corresponding 
to the {d— 1) Fourier modes in the range of the curl operator f C curl /, lying in (C/j^|)^. Thus, 
A{£^, A) = A-'^^''^(^, \) = is exactly the linearized dispersion relation for the constrained averaged 
system il.5\) . il.lCH) relevant to time- asymptotic behavior. 

Theorem ^ niay be regarded as partial justification of the WKB expansion. Roughly speaking, it 
states that if perturbed periodic waves exhibit coherent behavior near the unperturbed wave u, 
then this behavior is well-described by the constrained averaged equations ()1.5|) . (jl.lUj) . In one 
dimension, additional results of |OZ.2j case give rigorous sense to this statement in the form of 
detailed pointwise linear bounds under the assumption of spectral stability of the linearized operator 
about the wave. To establish a comparable long time result on behavior in multidimensions would 
be a very interesting direction for future investigation. Moreover, the new description of the solution 
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given in jSe.lj by modulation expansion might give a sufficiently good nonlinear Ansatz to carry 
out a complete nonlinear analysis, which was not done even in the one-dimensional case. Thus, it 
would be interesting to revisit also the one-dimensional setting of |UZ.2j from this new perspective, 
in particular, to resolve certain puzzling issues in the general, non-quasi-Hamiltonian case. 

An equally interesting direction for future investigation would be to rigorously validate the WKB 
expansion of Section |1] for the closely related small viscosity problem 



e — > 0, similarly as in 

iniwza-inMWTa for the viscous shock case. We note in this regard 

that a key ingredient in the |GMWZ.2j - |GMWZ.3j analysis, the conjugation lemma of |MZj assert- 
ing the existence of a coordinate change converting asymptotically constant-coefficient resolvent 
ODE to constant-coefficient form, has a straightforward analog in the standard Floquet construc- 
tion converting periodic- to constant- coefficient ODE; indeed, this construction is more natural in 
the periodic case. Thus, the whole Kreiss symmetrizer construction of |GMWZ?2] - |GMWZ.3j may 
be brought to bear also in the periodic case. Indeed, the situation is simpler since no shock front 
appears; as remarked in |OZ.2j . there is a much closer analogy to the constant-coefficient case. As- 
suming successful validation of the modulation equations, there is also the interesting question of 
secondary modulation, i.e., what happens after shock formation time for the averaged, hyperbolic 
equations? 

We point out further, as observed by Serre |Se.lj in the one-dimensional case, that the result 
of Theorem ^ is completely analogous to the corresponding relation established in |ZS| for the 
linearized dispersion relation associated with a perturbed viscous traveling front u = u{x ■ v — st), 
lim^^-l-oo = u-t, in which the WKB expansion corresponds to matched asymptotics joining an outer, 
hyperbolic solution and an inner, viscous profile. See |Z-kj . Section 1.3, for a formal derivation in 
the shock case, starting with the same rescaling — > (ex,ei); see also the related discussion of 
long-time vs. small- viscosity problems in |GMWZT| . Rigorous verification may be found in |Z-kj . 



As an immediate consequence of Theorem ^ we obtain the following two corollaries, yielding a 
necessary condition for low-frequency multi-dimensional spectral stability strengthening the one- 
dimensional version obtained in |OZ.lj . [Se.lj . 

Corollary 1.1. Assuming (H0)-(H2) and the nondegeneracy condition 



then for A,^ sufficiently small, the zero-set of D{-, ■), corresponding to spectra of L, consists o/n + 1 
characteristic surfaces: 





IC^MWZ.lj . 





iajiO + OiO, j = l,...,n + l, 



(1.13) 
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where aj(^) denote the eigenvalues of 



— Z^^J F)(M om ' ^^-^^^ 



d{M,QN) ' 



excluding {d—\) identically zero eigenvalues associated with modes not satisfying constraint hl.lU\) . 

Proof. Similarly as in as in the proof of the analogous Lemma 7.5 (ZSj in the shock wave case, 
assuming ()1.12|) . we may easily deduce ()1.13|) from (|1.9|) using Rouches Theorem. Defining 

pA5(A):=p-("+i)D(p|,pA), (1.15) 

for ^, A) € M X S"^^"*^ X C, we obtain a d-parameter family of analytic maps, converging as /? ^ to 
= Z\(^, •). Under assumption H1.12() . V^'^ = •) ^ A""^^ as |A| — > oo, hence, for p sufficiently 
small, VP'^ has n + 1 continuously varying roots A = aj{^,p). Defining aj(^) := |'?|aj(C/|Cl) we 
obtain the result. ■ 

Remark 1.1. Evidently, aj(^) are smooth in |^| for fixed ^, but in general have a conical singularity 
at ^ = when considered as a function of ^, i.e., daj/dS, is discontinuous at ^ = 0. 

Corollary 1.2. Assuming (H0)-(H2) and the nondegeneracy condition M.l^) . a necessary condi- 
tion for low-frequency spectral stability of u, defined as Re X < for D{^, A) = 0, ^ G W^, and 
1^, A| sufficiently small, is that the averaged system be "weakly hyperbolic" in the sense that it 
possesses a full set of real characteristics Xj{^) for each ^ € M*^, i.e., the eigenvalues of 

d{F^,Sf2ej) 



^ d{M,QN) 



are real. 



Remark 1.2. Condition (TTT^ . or equivalently (5/aA)"+^L>(0, 0) / 0, is a necessary condition 
for one-dimensional linearized stability |OZ.2j . while hyperbolicity is necessary for stability of the 
homogenized system linearized about a constant state. Thus, Corollaries 11.11 and 11.21 are analogous 
to results of ZS] in the shock wave case, stating that, given one-dimensional stability, stability of 
the inviscid equations linearized about an ideal shock is necessary for multi-dimensional stability 
of a viscous shock wave. 



Finally, we mention that, though the averaged system may in some cases be hyperbolic |OZ.lj . so 
far, only unstable periodic traveling-wave solutions have been found for viscous conservation laws. 
However, essentially only the single, 2x2 model of van der Waals gas dynamics with viscosity- 
capillarity in one dimension has so far been considered in detail |UZ.lj . |Se.lj - [Se.2j . and we see 
no obvious reason why a stable wave should not exist for other models. It would be extremely 
interesting to either find such an example, with the associated rich behavior described by the 
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modulation equations, or show that it can in no case exist. As suggested by Serre |Se.3j . a useful 
starting point might be to consider whether the averaged system ()1.5() might ever possess an entropy. 

Plan of the paper In Section [21 we recall the basic Evans function construction of [01 . In Section 
|31 we carry out the expansion ()1.7|) . and in Section |1J the multi-dimensional WKB expansion H1.5() . 
Finally, in Section [3 we carry out the proof of Theorem ^ by a calculation similar to the one used 
by Serre |Se.lj to treat the one-dimensional case. 

2 Preliminaries 

Without loss of generality taking S{u) = 0, N{u) = ei, u = u{xi) represents a stationary solution. 
Linearizing about u{-), we obtain 

vt = Lv := YiB^'^v^J,^ - YiA^v),^, (2.1) 

where coefficients 

B^^ := B^'^iu), A^v:=Dp{u)v-{DB^\u)v)u^, (2.2) 
are now periodic functions of xi. 

Taking the Fourier transform in the transverse coordinate x = (x2, ■ • • , x^), we obtain 

+ i{Y^B^^ikv)x,-iY,^'i3^- E B'^'^kijV. (2.3) 

where ^ = (^2, ■ ■ ■ ) ^d) is the transverse frequency vector. The Laplace transform in time t leads us 
to study the family of eigenvalue equations 

= (L^- - \)w = {B^^w'y - (A^w)' + i J2 B^^ijw' + i(^ B^^ikw)' 

-i^^A^^.w- B^%CjW-Xw, (2.4) 

associated with operators and frequency A G C, where ''' denotes d/dxi. Clearly, a necessary 
condition for stability of 1)1.1(1 is that (|2.4I) have no solutions w for £^ G M*^"^ and Re A > 0. For 
solutions of (|2.4|1 correspond to normal modes v{x,t) = e'^*e*^'^t(;(xi) of (|2.1|1 . 

The difficulty of our problem is due to accumulation at the origin of the essential spectrum of the 
linearized operator L about the wave as in the one dimensional case. Multidimensional stability 
concerns the behavior of the perturbation of the top eigenvalue, A = under small perturbations 
in ^. For study this stability, we use Floquet's theory and an Evans function |U] which not only 
depends on A but also on which corresponds to the phase shift and ^. To define the Evans 
function, we choose a basis {w^{xi,^, A), . . . , z/;^"(xi, ^, A)} of the kernel of L: — X, which is analytic 
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in A) and is real when A is real, for details see |UZ.llSeTl| . Now we can define the Evans function 
by 

u''(X,f,A) -e*^«iW(0,<f,A) 



(2.5) 



l<l<2n 

where .^i E M. Note that X^i is exactly 6 in |Se.lj . We remark that D is analytic everywhere, with 
associated analytic eigenfunction w'' for 1 < / < 2n. A point A is in the spectrum of if and only 

ifi)(A,o = o withe = (6,1)- 

Example 2.1. In the constant-coefficient case 

-i^A^^.w- Y B^\kijW-\w = Q, (2.6) 

an elementary computation yields 

D(A,0 = il'=i(e'''(^'«)^-e^«^^) 
where / = 1, . . . , 2n, denote the roots of the characteristic equation 

(^2^11 + ^1 + i Y B^^i, + i Y B^\k) 

-{iJ2A^^,+ Yl B^^iki,+\I))w = Q, (2.7) 

where w = q^^^^w. The zero set of D consists of all A and such that 

W(A,0 = iei(mod27rVX) 
for some I. Setting ^ = i^i in (|2.7j) . we obtain the dispersion relation 

det(-S« - iA^ - A/) = (2.8) 
where A^ = A^^, and B^ = ^^'efc^- 

Remark 2.1. If (A,|) = (0,0) then reduces to 

^{{B^^y^A^ - fx)w = 

giving n nonzero roots 

fi = Sj, W = tj, 

where Sj^tj are eigenvalues and eigenvectors for the matrix (i?^^)^^A^, and an n-fold root n = 0. 
Thus, -D(0, 0) = in the above example. We shall see later that this holds also in the general 
variable-coefficient case. 
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Remark 2.2. In the constant coefficient case, (|2.8j) yields expansions 

Xj{0 = 0-iaj{0 + Om, j = l,...,n, (2.9) 

for tlie n roots bifurcating from A(0) = 0, wliere aj denote the eigenvalues of A^. Thus we obtain 
the necessary stability condition of hyperbolicity, (7{A^) real. 



3 Evans function calculations 



Motivated by the example, we now find linearized dispersion relations for the variable-coefficient 
Evans function in the low-frequency limit. From now on, coordinatize v in the vicinity of ei by 



(5 = ((5^, . . . , (5'^) G R'^-'^. Note that differentiation of (EU) yields du = (0, 5). 

Lemma 3.1 r| Se.l p. Assumption (H2) is equivalent to 
(H2)' 



(3.1) 



span{ 



du 
ds 



du 
06 



,[w'],...,[w'%u'{0)} 



Proof. Immediate, using [wi] = 0; see |Se.lj for the one-dimensional case. 



3.1 Variational relations 

If = 0, (|2.4() equation becomes 

(Lo - X)w = {B^^w'y - (A^w)' -Xw = 0, (3.2) 

which is associated with the one dimensional stability problem studied in jUZ.llSeTTl . Recall that 
u is X-periodic in xi and the functions w^{xi,^, A), . . . , i(;^"(xi, ^, A) are in the basis of the kernel 
of — A. Following |OZ.lj . we normalize 

w^iO,i,X) = e„ iw^nO,i,X) = {B'^r'A^ef, 
^^"+^(0,e,A) = 0, (u;'^+^)'(0,|,A) = -{B^'r'e, 

for j = 1, . . . , n and all (^, A), giving in particular 

W(0, i, A) = 0, Lu;"+^ (0, i A) = ej for j = 1, . . . , n. (3.3) 

Plug the Taylor expansion of w{xi,(,, A) at the origin of (A,^) 

w'{;i A) = w'{; 0, 0) + Xw[{;0, 0) + ^ {; 0, 0)^, (3.4) 

+^(AMA(-,0,0)+2A5^t/;i^^,(.,0,0)e,+ Yl 0' 0)^.^,) + •• • 
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into (|2.4|) to find the identies: 

(Wy = 0, {Lwiy = w', {Lw{J = 2w{ (3.5) 

and 

{Lw\^J = (zA^wl + zA'^wl^ - tB^\wl)' - t{B'^w\J - tB'\w\^)' - i{B'^w\^y + 2i?^V), 
where Lw = B^^w' — A^w, and also: 

{LwlJ = i{A^w^ - B^\w^y - (B^^w^Y), i/1 (3.6) 

and 

{LwlJ = i{f{u)-B^\u)n' -B'^{u)ny, j^l (3.7) 
by using the definition of A^ in (|2.2|1 . In the Laplacian case B^^ = 5^, the latter identity simplifies 



to 



{Lwly = zP{uy. (3.8) 



Note that the functions u;^(xi,0,0), . . . ,w'^^{xi,0,0) are in the basis of the kernel of Lq. We omit 
(•,0,0) hereabove and denote [w] = w(X) — w{0). We also have = u', Lw^ = 0, [w^] = 0, 
w^dxi = 0, and moreover, 

[W]=0, [Lw{]= w^dxi, (3.9) 



and 



[Lwl] = 0, [Lwl] = 0, (3.10) 

Jo 

= £ {^A^^l - iB^\w\y - i{B^^w\y + wl^dx,, 



-iB''\wly - i{B^''wly + 2B^''w^)dxi. 



(3.11) 



In the Laplacian case B^'^ = dj., the last two identities simplify considerably, to 



[Lwl^,] = J^ {iA^wl + wl)dxi, 
[H€J = / {iA^wl+iA'wl)dx,. 



(3.12) 
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3.2 Connection to traveling-wave variations 



From (|3.3|) . we find easily that 

{-,0,0) = du/daj\u, w'^'^^ {-,0,0) = du/dqj\u for j = 1, . . . , n. 



(3.13) 



For example, taking the variation of traveling wave equation ()1.4() with respect to qj, we find that 
z = du/dqj satisfies 



B^^z' - A^z = -ej 



with z{0) = 0, so that Lqz = and z'{0) = —{B^^) ^ej as claimed. 
Further (see [071] . ^el\ ). 

2n 

wl = -du/ds + ^r^-^ 



a w 



n+l 



for a e M^", since Lo(—(?n/9s) = Low\{-,0,0) =u'andw\{0,0,0) = {du/ds){0) = 0, and, similarly, 
for j + 1, using Lqu;^^, (•, 0, 0) = Lo(i9n/9(5-')(-, 0, 0) = ip{u)' and it;^^, (0, 0, 0) = duldb^i^i) = 0, 



2n 



(•, 0, 0) = idu/d6^ + P\w^. 



n+l 



Alternatively, 



for 



^d1(.,0,0) = -du/ds, wli-,0,0) =idu/d6^ 



■w := w 



2n 2n 
n+l £=n+l j 



(3.14) 



(3.15) 



with 'w^{0,0) still equal to u'. We hereafter substitute for everywhere it appears, as we are 
free to do. (Recall, can be an arbitrary basis of the kernel of L.) 



3.3 Reduction of the leading part 



We rewrite the Evans function (|2.5|) as 



K(|,A) 



+ (l_e*^«i)V(0,e,A) 
+ (l-e^^«O(^«0'(O,e, A) 



(3.16) 



l<K2n 



and then multiply the second row in 1)3. 16(1 by B^^ and then subtract times the first one 



{detB'')D{X,^i,0 :-- 



K(C~, A)] + (l-e*^€i)^^'(0,e, A) 
Lw^iiX)] +(l-e*^«i)(Lti;')(0,e~,A) 



(3.17) 



l<K2n 
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At this point, we restrict for readibility to the Laplacian case B^^ = 6j.. The general case goes 
similarly. Then the Evans function D(X,(^i,0) becomes 



Iq det 



c(e, A) [w'] ...\ 

+ 0(|A|"+2 + [^^["+2) (3.18) 



with a nonzero number Fq (for details, see |OZ.llSeTT] ). where 

c(C, A) = X[wl] + J^^ KJ - iXCiu'iO) 



-X[du/ds] + ^iCj[du/d6j] - iXCiu'iO) (3.19) 



(3.20) 



is a homogeneous degree one polynomial, 

-iX^^XiLwlm - iX^,^^j{Lwl){0) 

= -X^ idu/ds)dxi+iXy2^j {A^{-du/ds) + {du/dP))dxi 
Jo Jo 

- E ^i^-t r A\du/d5^)dxi 

+iX^iX{Ldu/ds){Q) + X^i^ij{Ldu/d5^){Q) 

j-X j-X j-X 

= -AM {du/ds)dxi+iX^iA-{d/ds) f {u)dxi + {d/d6^) udx 
Jo ^ ^ ^0 ^0 



+iX^iX{Ldu/ds){0) + X(i^Cj{Ldu/d6^){0) (3.21) 



is a homogeneous degree two polynomial, and 



C^(e, A) := X[Lw{] + - iXCiLw\0) (3.22) 



are homogeneous degree one polynomials given by 
rX rX 



{d/daMxf u{x,)dx,+y2^j [ if{u{xi))dxi-iXUu'{0)-f{um) (3.23) 
^ Jo JO ^ 
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for £ = j = 2, . . . ,n, and 

id/dq,)(x r uixi)dxi + y2Cj [''if{uixi))dxi-iXCiiu'iO)-f\uiO))) (3.24) 
^ -'0 Jo ^ 

for I = n + j = 71 + 1,..., 2n. 

Thus, the leading order part of D near (^, A) = (0, 0) is the homogeneous degree (re + 1) polynomial 

' c{iA) ... \ 



Z\i(e,A) :=rodet 



c{i,\) C\i,\)... 



(3.25) 



with c, C, defined as above. In particular, the Evans function has a zero of order n + 1 at 

(A,e) = (o,o). 

4 Slow modulation approximation 

Next, we carry out a multi-dimensional version of the slow modulation (WKB) expansion in |Se.lj . 
Rescale (x, t) i— > (ex, et) in to obtain 

i j, k 



Let 

u\x,t)=u^ (x,t, '^^"^'^^ 



)+-^(x,t,^)+..., (4.2) 



e 

where y i— »■ vP{x, t, y) is a periodic function with dx(t> / 0. We plug (|4.2|) into (|4.1j) and consider the 
equations obtained by equating coefficients at successive powers of e. 

At order e~^, we have 

- sdyu^ + ^ujUjdy{P {u^)) - dyC^uo'^UjUkB^''{u^)dyU^) = 0, 

j j,k 

with 

'-■=-W¥y ^'-'^ 

which may be recognized as the traveling-profile equation after rescaling y ujy. That is, u^{y) = 
u{ujy) for a periodic profile of period X = uj^^, hence is periodic of period one, as described 
in |Se.lj . The quantities uj{x,t), s{x,t), h'{x,t) are the local frequency, speed, and direction of the 
modulated wave. 

At order e*^, we have 

d d 
dtu^ + dx, {f{u') - Yl B^\u'')uvkdyu^) = dy{. . .). 

j=l k=l 
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Taking the average with repect to y, and rescahng with y := ujy, we obtain 



dtMiu'>) + Y,d.,F\u^) = (4.4) 
i=i 



where 

fX 



= ^ / [P{u^) -Y,B'Hu>'^kdyu')dy (4.5) 

k=i 

is the averaged flux along orbit (now rescaled to actual period X), with 

^i/,F^ = (5M + Q)(nO), 
i 

by the profile equation. In the Laplacian case B^^ = S^., 1)4. 5(1 simplifies to 

F'{u^) = ^ £if'(^') - '^jin'')dyu')dy. (4.6) 

We have an additional d equations 

dt{nN){u^) + d^{nS{u^)) = (4.7) 

from the Schwarz identity dtdx(l) = dxdt(j), where d is the dimension of the spatial variable x. (Note: 
{il,N) may be regarded as polar coordinates for f]N.) 

Combining, we obtain finally the closed homogenized system 

dt{M, f2N) + dx, {F^, ^Sej) (4.8) 
j 

of the introduction, consisting of n+d equations in n+d unknowns. As discussed in the introduction, 
this should be supplemented with the constraint 

curl {QN) = (4.9) 

coming from the relation QN = Vx4>- 



5 Proof of the main theorem 

We now carry out the proof of Theorem^ restricting for readibility to the Laplacian case B^^ = S^.. 
The general case follows similarly. We want to see that the leading order part Ai of D, defined in 
(|3.25|) . is given by a (nonzero) constant multiple of A^"*^ times 

Z^(e,A) =det (a ^ \ u) + Y,^u:Xi, {u) 

where u denotes the orbit class of u, with i/ = (1, 0, . . . , 0) = N{u). 




14 



Recall the assumption (H2)' and p.lj) . The tangent space to P at u is the /3-projection of the kernel 
of 



Z(/3,5,7) :=/3° 



du 
ds 



2n d 

i=2 



1 



du 
W 



+ 7n'(0), 



and the tangent space to 7-" at u is the /3-projection of the kernel of 

2n d 



z(/3^/3^...,/3^^5^...,5^7) --=0 



du 
ds 



2 j=2 



du 
W 



+ ju{0). 



We relabel for f3^ since we will not use hereafter. Thus, 

2n d 



du 
9J 



2 i=2 



+ 7n'(0), /3gC^", (JeC^-Se 



(5.2) 



We easily compute (see 1)3. 3() for dQ) the differentials 

• (/3, (5, 7) = 7 = ax, dS-{/3,6,-f)=P^ =dS, dN ■ {13,6,-^) = {0,6f = dN, 

2n 

dQ ■ iP, 5, 7) = (/3"+\ . . . , P'n^ = -Yl z^'^^' = • ^) = -'^'^ = 



(5.3) 



where, following |Se.lj . we use the notation dG to indicate the extension to C^""'"'^ of a differential 
dG defined on the kernel of Z. (Note: this includes the extension from real to complex values, of 
which we shall later make important use in parametrizations ()5.14|) and H5.15() .') 

Likewise, XM = u{y)dy gives 
diXM)-iP,6,^)=jm+P' Y/y + Y.(^J^ w'dy + Y,^^ —dy = diXM), (5.4) 



thus determining 

dM = uj{d{XM) - MdX). 

Finally, variations 

d{XF^){u) = d I {f{u) - Uj{u)dyu)dy = d{XF^){u) 
Jo 

for J 7^ 1 have the simple form 



d{XF^){u) = dXP{u){0) + 
by 1)4. 6|) and periodicity of u, hence 



X 



^{u)dy, 



(5.5) 



(5.6) 



(5.7) 



dF^{u) =uj{d{XF^) - {dX)F^{u)) 

= ujjf{u){0)+uj [ dfiu)dy-uj-fF^u). 



(51 
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To find the variation for F^, note tliat, by the first-order traveling wave system H1.4() . 

Y^UjF^ = MS + Q, 

j 

so that + J2j ^jidF^) = d{MS + Q), hence, for z/ = (1, 0, . . . , 0), s = 0, 

dF^ {u) = d{MS + Q)-Y, ^^^^ (^) = ^^^^ + dQ-Y^ 5^F^ (u). (5.9) 

We may now compute the determinant ()5.1() . i.e., the determinant of the restriction to ker Z of the 
linear map 

H (3S^) = ( "^A-=( ^^^^^ + ^^^^ + iC,P'^e, \ 

A,€IP, ,1) ^^2^^(^^^^^)J • [xdM + iCoXC^{MdS + dQ -j:^^,6^Ff) + j:^_,,iLoXC,dF^ J 



which can be evaluated using an ingenious trick of |Se.lj as 
det(iJA,^lker2) = 



=CiC2(e,A)det ^'M IkerHHA,^' 



V 2 



(5.10) 



where Ci = det ^-^^Ikerz-i-) well-defined thanks to full rank of assumption (H2)', and 
independent of (^, A), by the corresponding property of Z, but C2(C, A) = detff^|]jgj.^j:^i^x and 



(A,0 



both depend on the specific dependence on (^, A) of the basis chosen for ker//^. 

Note that determinant 1)5. If) is in the first place defined only up to a constant factor depending on 
the parametrization of V, so that we need only take care of the (^, A) dependence of C2. Rewriting 

d d 
^i,€(A<5,7) = XiJ^OP -uj^jei) + Y,i^jP'^ej = (5.11) 
2 1 

as 

i^i^P'^ = XO^'J, (5.12) 

and 

Xlu5^ = -i^jOP^, j ^ 1, (5.13) 

and setting as in |Se.lj ^ 

P' = -Xp, -i = -iXiiP, (5.14) 



^ We make the inessential change p ~p for convenience in later calculations. 
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p E C, giving also the (new, multi-dimensional) relations 



5^ = iijp, 



(5.15) 



and leaving (3^ free for ^ 7^ 1, determines a choice of basis for keiH^, for which C2(^, A) has the 
simple form A'^"^. 

This fact is most easily verified by right-multiplying 

(Hi 



by the determinant one matrix 




/ a Xp \ 

02n-l 02n-l ^2n-l hn-l 

\ 15 iXilp ) 



aiX^ip — j3Xp = 1, to obtain 




N- 



1 



1 •= 



uj\Id-i ) 
)perator 

through the compositions H5.14() and H5.15() . and thus 

^(^,A) = CiA'^-MetiVs- 



where is the 2n x 2n matrix corresponding to linear operator | ^ 1 operating on (/), /3^, . . . , /J^") 



(5.16) 



Alternatively, we may observe that is full rank whenever A ^ 0. Observing also a posteriori that 
det is homogeneous degree n + 1, we may conclude that C(^, A) = det H^\-^qj.^^i^± as the ratio 
of n + d and n + 1 degree homogeneous polynomials must be a constant times A'^"^. This discussion 
repairs a minor omission in |Se.lj . where the dependence of C2 on (^, A) is not explicitly discussed. 

It remains to compute, under the compositions (|5.14l) . ()5.15|) . the 2n x 2n determinant detA^2) 
which, transposing first and second n-row blocks, may be expressed as 



det N2 = det 



XdM + iCdX^i{MdS + dQ- Y.j+x ^^F^) + Ej^i i^X^jdF^ 



(5.17) 



Substituting from H5.14() - (|5.15() and the variational formulae ()5.4|) - ()5.9|) . and expressing A^2 as a 
matrix taking (p, . . . , /3^") C^", we obtain, similarly as in |Se.lj . that detA^2 = ^A^"'^ is 
a)" times the determinant Z\i defined in 1)3. 25(1 . giving the desired relation Ai = oj^"^ X^^'^ A, and 
completing the proof. 
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Namely, the first line of (|5.17j) becomes 



2n 



-\p[du/ds\ + (3\w^] - iXiipu'{Q) + ^ i^jp[du/d6^ 



The second line of H5.17() becomes 



/ fX a ^ fX " _ \ 

+ V i^X^j f CJ7/^' (u) (0) + w / a/^' (u)dxi - O-fF^ (u) ) 

.yi \ Jo J 

( - du 2n X d du - - \ 

(2n d \ 

-MXp - Y P^Lw^ - Yl ^ijpP^{^) ) 

+ Y ^^^^i (^-OiXCipP{u){0) + w ^ df\u)dxi + u)iXiipF\u)^ 

\iXiim + >?j^ Q^dxi-XiY^jl ^dxi-X6^6/^(tl)(0)j 

2n / „x _ \ d. „x 

+ lD^/3M Ay w^dxi-iXCiLw^j+OiY^jJ df{u)dxi 



-Op 



/ ^ p 2n X <i j-X 



with H5.5p and other identities. Denoting by M the matrix in H3.25() for which Ai = det A/", we find, 
comparing term by term, that the first n rows of are equal to the first n rows of AA, while the 
last n rows of are equal to (D times the last n rows of TV. Thus, det = w^detA/" = uj'^Ai as 
claimed, and we are done. 
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